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Abstract. We report new simulations of cooling of compact stars containing quark cores and updated fits to the Cas 
A fast cooling data. Our model is built on the assumption that the transient behaviour of the star in Cas A is due to 
a phase transition within the dense QCD matter in the core of the star. Specifically, the fast cooling is attributed to 
an enhancement in the neutrino emission triggered by a transition from a fully gapped, two-flavor, red-green color¬ 
superconducting quark condensate to a superconducting crystalline or an alternative gapless, color-superconducting 
phase. The blue colored condensate is modeled as a Bardeen-Cooper-Schrielfer (BCS)-type color superconductor with 
spin-one pairing order parameter. We study the sensitivity of the fits to the phase transition temperature, the pairing gap 
of blue quarks and the time-scale characterizing the phase transition (the latter modelled in terms of a width parameter). 
Relative variations in these parameter around their best fit values larger than 10“'^ spoil the fit to the data. We confirm 
the previous finding that the cooling curves show significant variations as a function of compact star mass, which allows 
one to account for dispersion in the data on the surface temperatures of thermally emitting neutron stars. 

PACS. 97.60.Jd Neutron stars - 26.60.-c Nuclear matter aspects of neutron stars - 95.30.Cq Elementary particle 
processes - 74.25.Dw Superconductivity phase diagrams 


1 Introduction 

The central densities of compact stars (hereafter CSs) can be by 
a factor of few up to ten times larger than the saturation density 
of nuclear matter. In this range of densities the mean interpar¬ 
ticle distance may become of the order of characteristic size of 
a baryon, therefore neutrons, protons, and heavier baryons will 
eventually lose their identity, as their wave-functions start to 
overlap. Although the detailes of the mechanism of deconfine¬ 
ment are not well understood so far, it is plain that at sufficiently 
large density the baryonic matter will undergo a deconfinement 
phase transition to quark matter. Not only the transition itself, 
but also the phase structure of quark matter is difficult to ac¬ 
cess, because the densities in CSs will not be large enough for 
the perturbative QCD to be valid. Astrophysics of CSs offers 
an avenue to explore and constrain the possible deconfinement 
phase transition and the properties of dense QCD matter via 
modeling of CSs containing quark cores Qll]. 

The focus of this work is the effects of the deconfinement 
phase transition and the phase structure of dense QCD on the 
cooling of CSs. This discussion extends the previous work on 
the cooling of hybrid stars containing superconducting quark 
matter, including the interpretation of the rapid cooling of the 
CS in Cassiopea A as a phase transition within the QCD phase 
diagram ||^0l- Specifically, we report the results of new simu¬ 
lations of cooling of CSs and fits to the Cas A updated data a, 
which covers the 10 year period from 2003 to 2013. These 
data indicates an unprecedented fast cooling of the neutron star 
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in Cas A, which requires fast transient modeling of this ob¬ 
ject (as opposed to the familiar long-time-scale modeling of 
neutron stars, where characteristic timescales of variations are 
t > 100 y). The best estimate of Ref. a indicates a decline in 
the temperature of the star 2.9 + 0.9% over 10 years of obser¬ 
vation, which requires extremely fast transient cooling of the 
star over this period. The currently available data cannot be in¬ 
terpreted unambiguously; one reason is the bright and varying 
supernova remnant background, which makes a definitive inter¬ 
pretation difficult a . Another uncertainty arises from adopted 
constraints on the temperature fitting parameters and the uncer¬ 
tainties of the effective area calibration; for example, the appar¬ 
ent decline in the temperature can be compensated if one allows 
for variations in the emitting region size a^ In the following 
we leave aside the possible uncertainties in the interpretation 
of the data and focus on its theoretical modeling. 

A variety of models account theoretically for the Cas A 
data. Each class of models has its own specific trigger for on¬ 
set of rapid cooling around the age of the Cas A (f 300 
y). A class of nucleonic models (i.e., models of CSs contain¬ 
ing neutrons, protons, and electrons) attribute the rapid cool¬ 
ing to the onset of Cooper pair-breaking process in CS’s su¬ 
perfluid component [l^ . The rapid cooling 

occurs once the temperature drops below the critical temper¬ 
ature Tc^n - 10® K of phase transition of nucleons to the su- 
perluid state. This allows for additional neutrino emission via 
Cooper-pair breaking and formation (PBF) processes from S - 
wave |[l3l [13 Ell ull E^l and E-wave condensates ll^ 1^ . 
These minimalistic models do not require any additional new 
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physics beyond the standard scenario which is supplemented 
by pair-breaking processes jUt]. 

The Cas A data was also explained in Ref. within a 

model where modified Urea and bremsstrahlung processes are 
enhanced compared to rates used in the minimalistic models 
quoted above by several orders of magnitude due to a softening 
of pionic modes, as initially discussed in Ref. Rotational 
changes in star’s composition may induced Urea process as the 
star slows down, which in turn can cause a Cas A type fast 
transient 1^ . which however would require very fine tuning 
of parameters, as already discussed in Ref. JH]. 

Color-superconductivity was included in the simulations of 
CS to account for Cas A data in Refs. In Ref. the 

phase transition from purely nucleonic tc^urely quark matter 
occurs via mixed phases, whereas Ref. assumes sharp in¬ 
terface between the two; which scenario is realized depends on 
the poorly known surface energy between the nucleonic and 
quark phases. We explore here the second possibility. 

This paper is organized as follows. In sect.|2]we discuss the 
relevant features of the phase diagram of dense and cold QCD. 
Section |3] describes the physical input needed for cooling sim¬ 
ulations, including neutrino emissivities of various phases and 
specific heats. The results of numerical simulations are pre¬ 
sented in sect. 0] together with the fits to the Cas A data. A 
summary and conclusions are given in sect. |5] 

2 The phase structure of dense QCD and 
massive compact stars 

This section provides the background information on the QCD 
phase diagram, which is necessary for the understanding of the 
model of cooling of CSs presented in the following sections. 
Readers interested in the astrophysical aspects of the modeling 
can proceed to sect. [3j those interested in the results can go 
directly to sect. 0] 

We consider a class of hybrid CSs , which contain dense 
quark cores surrounded by a nucleonic envelope, with a sharp 
phase-transition interface between these phases. Stars with such 
hybrid structure naturally correspond to the most massive mem¬ 
bers of the sequence of stellar equilibria modelled with an equa¬ 
tion of state which contains a phase transition from nucleonic to 
quark matter. These massive members must be heavy enough to 
account for the recently observed two-solar mass pulsars PSR 
J1614-2230 and PSR J0348-1-0432 103[H. 

If quarks in the deconfined phase of the star are unpaired 
(i.e. do not form Cooper pairs) they will cool the star by neu¬ 
trino emission rapidly to t emp eratures well below those ob¬ 
served in Cas A SUllllSlllllMl- However, cold quark mat¬ 
ter is expected to be in one of the conjectured superconduct¬ 
ing phases due to the attractive component in the gluon ex¬ 
change quark-quark interaction |03- The y6-equilibrium and 
strange quark mass shift the Fermi surfaces of up and down 
quarks apart. As a consequence some non-BCS phases can now 
emerge. Such non-BCS phases include the gapless two-flavor 
phases 03 OF for example, the crystalline color-super- 
conducting phase 103 03 ■ The crystalline-color-superconduc¬ 
tivity itself has a multitude of realizations, which differ by the 
way the translational symmetry is broken by the condensate of 
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Fig. 1. Phase diagram of color-superconducting quark matter in the 
temperature (T) and isospin imbalance a = (nj - plane, 

where and «„ are the number densities of d and u flavors of quarks. 
The arrows schematically show the path in the phase diagram during 
the cooling of a star. 


Cooper pairs carrying finite momentum. For concrete calcula¬ 
tions we will assume below the so-called Fulde-Ferrell phase 
(hereafter FF phase), which is simple to model, but is gen¬ 
eral enough to preserve a key feature of the crystalline phases, 
which is the existence of gaple ss modes on the Fermi surfaces 
of up and down quarks 108103.1^1^ . 

Our conjecture concerning the behaviour observed in Cas 
A is based on a robust feature of the phase diagram of the 
two-flavor color superconducting matter, (and, in fact, any two- 
component fermionic systems with attractive interaction). Just 
below the critical temperature and at not too high isospin asym¬ 
metries, measured by parameter a - (iid - nu)f{nu + rid), where 
lid and are the number densities of d and u flavors of quarks, 
the rotational/translational symmetries are unbroken, i.e., the 
fully gapped 2SC phase is favored (see Fig. [TJ. At lower tem¬ 
peratures a phase transition to a less symmetrical phase (such 
as the FF phase) becomes favorable at larger asymmetries, but 
not large enough to destroy the superconductivity completely. 
(This would correspond in Fig. [T]to a > 0.2). Note that the 
phase diagram shown in FiglT] is valid only in the weak cou¬ 
pling limit, which is the case in the bulk of CSs . (If the cou¬ 
pling is large enough the ensemble makes a transition from 2SC 
to Bose-Einstein condensed phase without entering crystalline 
phase; see, for example. Refs. 143. 03l 03 for a generic and 
most complete phase diagram of a fermionic system with im¬ 
balance). 

We next explore the consequences of the transition from the 
symmetrical two-flavor BCS phase to the crystalline phase (in 
our case the FF phase) for the cooling of CSs and the CS in 
Cas A. We anticipate such a transition in CSs because a newly 
born proto-CSs have temperatures of the order of several tens 
of MeV and matter in nearly isospin symmetrical (a ^ 0) state. 
Within short period of time (minutes to hours) the CS cools 
to temperatures of the order of 1 MeV and becomes strongly 
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isospin asymmetrical as j6-equilibrium among quarks and elec¬ 
trons is established. The arrangement of phases shown in fig. [T] 
in the case of j8-equilibrated matter is shown in fig.|2] which is 
adapted from Ref. lliHl . A shell of quark matter phase in the star 


Table 1. The table lists the central density Pc.w, quark core radius Rq, 
quark plus hadronic core radius Rq+h, the isothermal core radius R„, 
the star radius R, and the masses of the quark core Mg, the hadronic 
core Mfj, the crust M^.^, and the total mass M. The density is in units of 
10*^' g cm“^, the radii in kilometers, and the masses in units of Mq 0]. 



Dg [fm 


Fig. 2. Phase diagram of the same phases as in fig. [T] but in the 
temperature and baryon density (ng) plane for quark matter in /?- 
equilibrium with electrons. The hadronic phase is shown schemati¬ 
cally on the left. 


with fixed density will traverse the T-a plane from upper left 
corner to the lower right corner as shown by arrows in fig. [T] 
The same path is shown in fig.l^in the temperature and density 
plane. 

The transition temperature from the 2SC to the FF phase is 
denoted by T* and will be treated below as a free parameter. 
Our fits to Cas A data below will demonstrate that the values 
obtained from theoretical models are larger than those obtained 
from fits. This uncertainty is acceptable as the magnitude of 
the pairing interaction among quarks in the non-perturbative 
regime of QCD is model dependent and generally not well- 
know. 

Finally, assuming that the cross-flavor pairing occurs be¬ 
tween the green and red colored quarks we need to address the 
pairing among the remaining blue colored quarks. The strength 
and the flavor content of pairing among blue quarks is model 
dependent Obviously blue-up and blue-down pair¬ 

ing is suppressed by a large mismatch between the Fermi sur¬ 
faces of the up and down quarks. Pairing is more likely in the 
color 6s and flavor 3s channel, which is same-flavor and same- 
color pairing, therefore it is not affected by the flavor asymme¬ 
try 11431] . Calculations show that the gaps are in the range 10-100 
keV. Again because of uncertainties envolved in the interaction 
generating the BCS state of blue quarks we treat the blue quark 
gap zife as a parameter of our model. 
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3 Physical input into simulations 

3.1 Equilibrium structure of hybrid stars 

Prior to computing the time evolution of CS models with color 
superconducting phases discussed in the previous section we 
need to construct the equilibrium models of such stars. Two 
classes of models were developed and used by us in the studies 
of color superconducting quark matter. Both classes are based 
on equations of state (hereafter EoS) with a sharp phase transi¬ 
tion between the nuclear core and quark phases; the first class 
features crystalline color superconducting quark matter in the 
quark core and nuclear envelope, which is described by a rel¬ 
ativistic density functional model with a stiff parameterization 
of the EoS The second class of models contains 2SC 

phase discussed in the previous section and the color flavor 
locked (CEL) phase at high densities ifsollsil] . Both models pro¬ 
duce stable sequences of color superconducting CSs with max¬ 
imal masses above the lower observational bound mentioned 
above. We model the high-density matter EoS within the effec¬ 
tive Nambu-Jona- Lasinio four-fermion interaction model sup¬ 
plemented by pairing interaction between quarks. The nuclear 
and quark equations of state were matched by a Maxwell con¬ 
struction. If one uses alternative constructions mixed phases 
will appear, which we do not discuss here (see Refs. 

We will assume that strange quarks because of their large mass 
are not important at the densities relevant to our models. The 
nucleonic phase itself can contain non-zero strangeness in form 
of hyperons. But the compatibility of hyperons with large masses 
of pulsars needs special attention. Indeed their appearance soft¬ 
ens the EoS and the maximal masses of neutron stars are lower 
than those observed. This “hyperonization puzzle” is a funda¬ 
mental problem on its ow n rig ht and will not be discussed here 
(see Refs. ll53ll5^l55ll56ll57Ll58ll5^ for recent discussions of 
this problem). Below we will exclude the hyperons from the 
consideration. 

The parameters of our models are listed in table [1] The 
models with low central densities (pc < 3po) have masses M < 
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1.85M0 and are purely hadronic. The models with larger cen¬ 
tral densities contain quark cores, with the size of the core in¬ 
creasing with the central density. 


3.2 Neutrino emissivities in quark phases 


We consider first unpaired quark matter. Because quarks are 
ultra-relativistic the Urea process works without any kinemati- 
cal restrictions (as opposed to baryonic matter). For two flavor 
quark matter consisting of u and d flavors of quarks the main 
modes of neutrino emission are given by 

d ^ u + e + V, u + e d + V, (1) 


where e stands for electron, v and v electron neutrino and an¬ 
tineutrino. The corresponding emissivity (energy output per unit 
volume and time) was calculated to leading order in the strong 
coupling constant as ; the emissivity per quark color is given 
by a, 

914 

ep = epdPuPeCUsT^, ( 2 ) 

where G is the weak coupling constant, 0 the Cabibbo angle, 
and pd, Pu, and pe are the Fermi momenta of down quarks, up 
quarks, and electrons. 

The quark pairing modifies the temperature dependence of 
process (|2]i. If the quark condensate was a BCS-type super¬ 
conductor, the emission would have been suppressed linearly 
for T ^ Tc and exponentially for T <sc Tc, where Tc is the 
critical temperature. However, in the gapless superconductors 
there is an additional new scale 6p = {pd - Pu)!'^, where pup 
are the chemical potentials of light quarks. Depending on the 
relative ratio of the two scales in the problem, namely 6p and 
Aq defined as the pairing gap for dp - 0 the suppression of 
emissivity is qualitatively different. In terms of the parameter 
^ = Aoidp one finds that in the ^ > 1 regime the Fermi sur¬ 
faces of quarks are gapped and the emissivity is suppressed 
as in the case of BCS superconductors iH. In the opposite 
case f < 1 the Fermi surfaces have nodes and particles can 
be excited around these nodes without energy cost needed to 
overcome the gap. In the case of the FF phase the shift in the 
chemical potential is replaced by a more general function - the 
anti-symmetric in the flavor part of the single particle spec¬ 
trum of up and down quarks. As a result the parameter dp is 
replaced by {ed{Q) - e«(6)]/2, where Q is the total momentum 
of a Cooper pair and eujd are the spectra of the u and d quarks. 

It is convenient to use below the generic parameterization 
of the suppression factor of the quark Urea process given by 


<T,)^2f{0ep, 


no- 


1 

exp[(^-l)f-i] + r 


(3) 


A time-independent constant ( excludes the possibility of the 
phase transition between the above mentioned phases. There¬ 
fore, following Ref. il , we adopt temperature dependent (and, 
therefore, time-dependent) parameter f(7’). We use the follow¬ 
ing parameterization 


(4) 


where is the initial (pre-transient in the Cas A case) value, 4^ 
the constant change in this function, and the function g{T) de¬ 
scribes the transition from the initial value to the asymptotic 
final value - f,- -A(. The transition is conveniently modeled 
by the following function 


giT)^ 


1 



(5) 


which allows us to control the temperature of transition, con¬ 
trolled by T*, and the smoothness of the transition, controlled 
by the width w. 

For the blue quarks, which are not involved in the FF pair¬ 
ing, we assume that these are paired as in the BCS phase, as al¬ 
ready discussed in the previous section. In that case their emis¬ 
sivity in the superconducting phase is very approximately re¬ 
lated to the emissivity in the normal phase 


4iT <T,)^^p{T>Tc)Jb 



nAb \ 


( 6 ) 


where the function Jb is some power-law function and n is an 
integer of order of unity. We assume 7^, = 1 as it is of order 
unity when T ~ Ab and is dominated by the exponential factor 
for T <s: Ab- The pairing gap Ab will be treated as a fit parameter 
below. 


3.3 Neutrino emissivities in hadronic phases 

Our treatment of the hadronic phases follows the standard pic¬ 
ture of slow cooling, which is supplemented by the pairing¬ 
breaking processes in the superfluid phases. Because of low 
proton fraction (xp ~ 0.05) in our models the dominant neu¬ 
trino emission channel in the unpaired phases of hadrons is the 
modified Urea process. Additional smaller contribution from 
the modified bremsstrahlung process is also included. We use 
emissivities derived for the free pion-exchange model of strong 
interaction 1^ . but reduce them by a constant factor 5 to ac¬ 
count for short-range repulsive component of the nuclear force E 
M . The effective mass corrections are included as well. 

Below the critical temperatures (Tc) for neutrons and pro¬ 
tons the emissivities of these processes are suppressed, at asymp¬ 
totically low temperatures exponentially. In general the emis¬ 
sivity of the modified Urea process in the superfluid phase can 
be expressed via its rate in the normal matter as follows 


epmodiT < Tc) - epmodiT > Tc) 

mA„(T) + kAp{T) 


^ Jpmod eXp I 


(7) 


where m and k are integers and 7^ mod is a power-law function 
of ApjT and AnjT, where An,p are the pairing gaps of protons 
and neutrons. The emissivities of the nucleonic superconduct¬ 
ing phases were derived in Ref s. ll^l64r for both modified and 
direct Urea as well as in Ref. 11^ 1^ for the direct Urea pro¬ 
cess using different methods. Below, for sake of simplicity, we 
assume m - k - 1 as well as 7^ mod = 1- This approxima¬ 
tion captures approximately the low-temperature asymptotic 


an^^i-A^giT), 
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reduction of the emissivity, but is a rather crude approxima¬ 
tion at r < Tc- More detailed expressions are not needed for 
the description of cooling of hybrid stars, which is dominated 
by the processes involving quarks. The neutrino emission pro¬ 
cesses from Cooper PBF start to contribute below respective 
Tc Sill. The neutral vector current processes are strongly 
suppressed by multiloop contribution to the response function 
of S -wave paired condensate The axial-vector 

emission can be taken at one-loop level to a good accuracy. For 
P-wave paired neutrons the rates are unaffecte d by the multi¬ 
loop processes and we use the results of Refs. ll2lLl^ . Given 
the zero-temperature values of the nucleonic gaps A{Q), their 
finite-temperature values are well approximated by the asymp¬ 
totic expressions given in Ref. (as well as fits to tables in 
that work 1^ 1 

A{T) f 1 - 0 < T < 0.5, 

zl(0) “ \ V3.016(1 - t) - 2.4(1 - t)2 0.5 < t < 1, 

where A(0) is the pairing gap in the zero-temperature limit, 
T - TjTc is the temperature in units of the critical tempera¬ 
ture and y - 1.781. The formulae (O reproduce the solution 
of the BCS gap equation with zero-range interaction within a 
percent accuracy. Crustal neutrino emission contributes to the 
cooling of the star in the final stages of neutrino cooling era via 
the electron neutrino bremsstrahlung emission on nuclei in the 
crusts ill. Ions in the crust may form a fluid or a solid. The 
abundances of impurities, if high, could lead to additional the 
neutrino emission. In the fluid or impurity dominated crust the 
emissivity scales as T® and we adopt this option for temper¬ 
atures relevant for that particular cooling phase {T < 10®/T). 
If crystalline lattice is formed at the temperatures of interest, 
then the emissivities would scale approximately as T® (classi¬ 
cal crystal) or as (quantum crystal) and the bremsstrahlung 
emissivity is parametrically suppressed iH. 


3.4 Specific heat 

The temperature dependence of the specific heat of normal Fer- 
mi-liquids (both relativistic and non-relativistic) is given by the 
linear law cy - aT, where the coefficient depends on the abun¬ 
dance of given species (which is the same as the Fermi momen¬ 
tum). We apply this formula to the leptonic component of the 
star and unpaired hadrons and quarks. 

As well known, in ordinary BCS superconductors the spe¬ 
cific heat as a function of temperature experiences a jump at 
the corresponding critical temperature and then decays, at low 
temperatures exponentially. For S -wave superconductors one 
can apply the a asymptotic expressions given in Ref. (as 
well as fits to tables in that work to model the ratio of the 
specific heats in the superconducting cj and normal cm phases 

( (127r/7)(2yT)-3/2e-'^/(rb Q < r < 0.3, 

= i -0.24422 -H 0.255292T (9) 

<^n(Tc) 1 h-2.43077t" 0.3<t<1, 


To find the specific heat of the red-green condensate we 
need an expression for the critical temperature as a function of 
the mismatch b/r. We use the following formula 


TciO - Tco ^1 - ( 10 ) 

where ^ (jud + Aq = A(^ = 0), Tco = Tc(^ = 0), and 

6(0 = (nd - n„)/(nd + riu). 

As in the case of emissivities, the availability of gapless 
fermions in the case f < 1 changes the BCS behaviour of the 
specific heat as well. To model this regime we apply, in analogy 
to emissivities, the following formula 


T < Tc) 




( 11 ) 


where is the specific heat of red-green unpaired quarks, 
taken as that for noninteracting quarks and is the specific 
heat of pair-correlated quarks. 


4 Results of numerical simulations 

4.1 General considerations 

The thermal evolution code employed in our study uses the 
isothermal core approximation which is valid for timescales 
that are larger that those which are required to dissolve tem¬ 
perature gradients by thermal conduction. In the hadronic core 
the thermal conductivity is dominated by the electron trans¬ 
port if baryons are superfluid iFt^I and baryons if these are non¬ 
superfluid iIt^ . The characteristic times-scales for thermal re¬ 
laxation are of the order of f ~ 100 y. In the 2SC phase likewise 
the transport is dominated by electrons and blue-quarks iH. 
the early stage of thermal evolution of 2SC phase and the time- 
scale of its thermal relaxation remains to be studied; here we 
assue that the entire (quark plus hadronic) core is thermally re¬ 
laxed at the start of simulation. Note that for reasonable initial 
temperatures the cooling tracks exit the non-isothermal phase 
and settle at a temperature predicted by the balance of the dom¬ 
inant neutrino emission mechanism and the specific heat of the 
core independent of the earlier evolution phase. 

The isothermal core is defined by the transition density pu = 
10'® g cm“^. At lower densities the envelope maintains temper¬ 
ature gradients throughout the complete evolution. The physics 
of thermal transport can be encoded in a relation between the 
interior and surface temperatures. The models of envelopes pre¬ 
dict the scaling T^ = gsh(T), where Ts is the surface tempera¬ 
ture T is the the isothermal core temperature, gs is the surface 
gravity, and h is a function which depends on T, the opacity of 
crustal material, and its EoS. This can be written as 117^ 


= (aTgy^gl/^ (12) 


where Tc is the critical temperature of phase transition. Equa¬ 
tion (|3 was applied uniformly to neutron, proton, and blue- 
quark condensates. See Ref. 11^ for possible effects of P-wave 
superfluidity of neutrons in the core on their specific heat. 


where Tg the isothermal core temperature in units of 10® K, 
Ts6 is the surface temperature in units of 10® K, a and p are 
constants that depend on the composition of the star’s atmo¬ 
sphere, gs is the surface gravity in units of 10'"' cm s"^. We 
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Fig. 3. Dependence of the surface temperature of a CSs given in K 
on time in year. The lines correspond to different phase transition tem¬ 
perature T* for fixed values of the width w = 0.2 and the blue quark 
pairing gap zlj = 0.15 (both in MeV). The labels correspond to values 
of T* in MeV. The Cas A data is shown by squares. 


use p - 0.55, which lies between the values for the purely- 
iron {p - 0.5) and the fully accreted {p - 0.61) envelope and 
a - 18.1, which is appropriate for an accreted envelope A 
Our fit results are sensitive to the p parameter, but not to the 
a parameter. With these ingredients at hand we have numeri¬ 
cally integrated the thermal evolution equation as described in 
Ref. @1. 

The numerical simulations were carried out for all models 
listed in Table [1] but we concentrate below mostly on the rep¬ 
resentative model with M/Mq = 1.93. Three parameters were 
varied; (a) the transition temperature T* from the 2SC to the 
crystalline phase; (b) the blue quark gap Ah', and (c) the width 
of the transition w. The remaining parameters which model the 
temperature (time) dependent part of the ( function were held 
fixed at values = 1.1 and - 0.2. The constant zero- 
temperature value of the red-green gap is fixed at the value 
Arg = 60 MeV. Figure [3] shows the dependence of the (red- 
shifted) surface temperature on time for M/M© = 1.93 stellar 
model and for different values of the T* (as indicated in the 
figure) with w and Ah fixed. In the limit T* —> 0 the quark 
core does not influence the cooling because the neutrino emis¬ 
sion from the red-green and blue condensates is suppressed for 
T <sc Arg and T Ah- For high values of T* (0.2 MeV in the 
figure) the transition to the FF phase occurs early in the evolu¬ 
tion of the star, therefore enhanced neutrino emission cools the 
star rapidly below the value observed for the CS in Cas A. It 
is further seen that by tuning the phase-transition temperature 
to the value T* - 0.009385 MeV the “drop” in the temper¬ 
ature of the model can be adjust to the observed temperature 
and age of the CS in Cas A. Thus, if the current interpreta¬ 
tion is correct, then the CS is a massive compact star (so it 
features a quark core) and it undergoes currently a phase tran¬ 
sition from the 2SC to the FF phase (or some other gapless 
phase). Note that the deduced phase transition temperature is 



Fig. 4. Same as in fig.[3]but for different values of the width w for fixed 
T* = 0.009385 and Ai, = 0.15 MeV. The labels correspond to values 
of w in MeV. 


of the order of 0.01 MeV and is small compared to the scale 
set by the Arg. This could be the feature of the equilibrium 
phase diagram of the superconducting quark matter or alter¬ 
natively may reflect the non-equilibrium aspect of the phase 
transition, i. e., long-lived metastable 2SC phase (in parallel to 
super-heating or super-cooling of superconducting states, know 
from the physics of ordinary superconductors). We conclude 
that with a proper choice of the value of the transition tempera¬ 
ture we can account for the observed “drop” in the temperature. 

Next we fix the parameters T* and Ah and vary the width of 
the transition w, which controls the “smoothness” of the tran¬ 
sition (fig.lH. For small values of w (0.1 in the figure) the tran¬ 
sient is steeper, as expected, but is also delayed (i. e. the drop is 
shifted away from the Cas A location to later times). Inversely, 
large w make the transient much smoother and also the start of 
the phase transition is shifted to the left, i. e., to earlier times. 
Thus, small (by a factor of 2) variation of the parameters w 
and T* can cause substantial shifts in the cooling curves, i. e., 
the overall fits to the Cas A data are sensitive to the values 
of these parameters. Finally we fix the parameters T* and w 
and vary the value of the gap of blue quarks Ah (fig. |5]l. In¬ 
creasing the magnitude of Ah shits the transient to later times, 
but also the pre-transient temperatures increase because of the 
exponential suppression factor in eq. (|6]l is less effective. In¬ 
versely, decreasing the value of Ah shifts the transient to earlier 
times and decreases the transient temperature. This behaviour 
is consistent with the fact that the cooling for heavy models 
is dominated by the red-green and blue quark components. In 
the 2SC phase the cooling is controlled by the Urea process on 
blue quarks, which is suppressed for late times and low temper¬ 
atures, see Eq. (I6]l. Unless the value of Ah is too small the blue 
quarks dominate throughout the neutrino cooling era. At later 
stages t ~ 10^ y the photon emission becomes the main cool¬ 
ing mechanism (this is indicated by the change of the slope of 
the cooling curve). During and after the transient, the gapless 
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Fig. 5. Same as in fig.[3]but for different values of the blue quark gap 
Ai, for fixed T’ - 0.009385 and w = 0.2 MeV. The labels correspond 
to values of /!;, in MeV. 



Fig. 6. Same as in fig. [3] for fixed T" = 0.009385 but for different 
masses of CS. The labels correspond to masses of the stars in units 
M/Mq. The references in the data are given in fig. 3 of Ref. |0|. 


phase emits neutrinos on red-green quarks; their luminosity is 
larger by factor two - the number of involved flavors - com¬ 
pared to blue quark luminosity. In this case again the transition 
is to the photon cooling at f ~ 10^ y. An important conclusion 
is that for massive CS the neutrino cooling is completely dom¬ 
inated by quarks and the details of the modeling of hadronic 
neutrino cooling are irrelevant. This implies, obviously, that the 
quality of the fits to the entire population of thermally emitting 
pulsars is not correlated to our fits to the Cas A data, as this 
population needs to include purely hadronic neutron stars with 
M < I.SMq, which will cool much slowly. 

We show in fig. the variation of cooling curves with the 
mass of the CS holding the “best” values of the parameters fit- 



t [yr] 


Fig. 7. Dependence of the surface in units of 10* K temperature on 
time for several temperatures of the phase transition T* at fixed w = 
0.2 and d* = 0.15 (both in MeV). The dots with error bars show the 
Cas A data taken from (5. 


ting Cas A fixed. It is seen that there is a substantial variation 
in the cooling speed of CSs with the mass; the heavy stars are 
cooling much faster than the light stars, for the neutrino emis- 
sivity of the hadronic phases is less effective that those of the 
quark phases. Note also the inversion of the temperatures in the 
photon cooling era; here the heavier M/M© = 2.05 star is hotter 
than the lighter M/M© = 1.93 star. The band lying between the 
cooling curves of extrem mass object, i.e., 1.1 < M/M© < 2.05 
in the T-t plane can be covered by changing the central den¬ 
sity (and therefore the mass) of the model and the experimental 
data can naturally be accommodated. (The few outliers, which 
are hot at late stages of evolution t ~ 10® y, may reflect differ¬ 
ent surface composition, than assumed in our simulations, e.g., 
composition containing hydrogen). Thus, we conclude that the 
spread in the temperatures of thermally emitting neutron stars 
can be attributed to the different masses of these objects; the 
lighter ones not containing quark matter cool more slowly than 
the heavier stars having quark cores. Of course, variations in 
other factors, such as the strength of the magnetic field, the 
surface composition, etc can also contribute to the dispersion. 

4.2 Fits to the Cas A data 

The Cas A data spanning the decade 2003-13 have been used 
for the fits. It is assumed that the star is already in the transient 
stage, as there is no evidence of the entry into this phase from 
a slower cooling stage, as well as no sign of the exit from the 
fast cooling mode. The parameters were adjusted to obtain the 
location of the Cas A data in the T-t diagram as well as the 
slope of the transient. The best fit is obtained with the parame¬ 
ter values T* - 0.009385, w = 0.2 wAAh - 0.15 (all quantities 
in MeV). The sensitivity of fits was tested against variation of 
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Fig. 8. Same as in|7] but for a range of values of w indicated in the plot 
and for fixed values of T" = 0.009385 and At, = 0.15. All parameter 
values are in MeV. 



t [yr] 

Fig. 9. Same as inj?] but for a range of values of A), indicated in the 
plot and for fixed values of T* = 0.009385 and w = 0.2. All parameter 
values are in MeV. 


one of them at fixed values of the other two. Figure|2]shows the 
effect of variations of T* on the transient behaviour. We quan¬ 
tify the sensitivity of the fits by defining the relative deviation 
as 6T*IT*, where 6T* is the range of values of T* compatible 
with the data; we find that ST* jT* ^ 10“^. Figure |8] shows the 
effect of variation of the width w at fixed values of other two 
parameters. In this case we find that dwfw = 1.5 x 10“^, i.e., 
the relative deviation is in the same range as for the transition 


temperature T*. Figure|9]shows the variation in the value of the 
blue quark gap zl;,; in this case we find that dzlfo//!;, ^ 3.3x10“^. 

Thus we conclude that the fitted values are quite sensitive 
to the precise values of parameters and small deviations at the 
level of 10“^ can spoil the fit to the Cas A data. Extrapolating 
from the present fits one would predict that the CS in Cas A 
would continue the fast cooling up to temperatures somewhat 
below 10^ K. Clearly, an observation of an exit from the tran¬ 
sient would provide further strong constraints on any model of 
the behaviour of the CS in Cas A and current fits should by 
updated to account for that effect. 

5 Summary and conclusions 

The observation of two-solar mass pulsars in binary systems is 
a strong evidence that the EoS of dense matter must be stiff and 
that the densities in the centers of compact stars can be large 
enough (several times the saturation density) so that the thresh¬ 
old density for transition to quark matter can be reached. This 
motivates the studies of the dense phases of quark matter un¬ 
der the conditions expected in compact stars (charge and color 
neutrality, near-equilibrium with respect to weak processes). 
We have concentrated here on CSs that have sharp interface 
between the hadronic and quark matter phases, in which case 
the massive members of the sequences (M > 1.85Mo) contain 
progressively larger amount of quark matter whereas the lighter 
members are purely hadronic. We have evolved these models in 
time and followed the changes in the core and surface temper¬ 
ature under assumption that the core is isothermal, which is a 
valid approximation, except for the first 100 y following the 
birth of the star. 

The generic phase diagram of imbalanced fermionic sys¬ 
tems, which also includes the flavor asymmetric pairing among 
quarks, implies that at high temperatures the pairing is in the 
BCS type phase, i.e., both Eermi surfaces of quarks are gapped 
despite of the mismatch between the Eermi surfaces of u and d 
flavors. At lower temperatures a transition from the BCS to a 
generic gapless phase must take place. This gapless phase will 
have some type of spatial modulation and we considered the 
EE type simple realization of such phase. The neutrino emis- 
sivity of such phase must be much larger than the emissivity of 
the 2SC phase because of the existence of gapless excitations. 
We have modelled the emissivities of both phase in terms of a 
simple parameterization, which however takes into account the 
fact that as the temperature is lowered there is a phase transition 
from the 2SC to EE phase. 

The rapid cooling of the CS in Cas A can be accounted 
for via the phase transition described above; we stress again 
that the phase transition takes place within the phase diagram 
of QCD and is the consequence of the ordering of various su¬ 
perconducting phases in the temperature, density, and isospin 
spaces. Such ordering was observed in numerous studies of 
other imbalance systems such as the hadronic superfluids or 
ultra-cold atomic condensates. 

The Cas A data can be fitted by varying the phase transi¬ 
tion temperature T* at fixed value of the remaining physical 
parameters - the width of the phase transition (or, equivalently, 
its duration) and the gap in the spectrum of BCS-paired blue 
quarks. It turns out that the relative accuracy of order 10“^ is 
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required to fit the data, i. e., larger deviation will spoilt the fit. 
We also tested the sensitivity of the other two parameters to the 
variations when the remaining two parameters are fixed. The 
relative accuracy is again in the range 10“^. Assuming that the 
present model of cooling of CS in Cas A is correct, we immedi¬ 
ately conclude that only the massive members of CS sequences 
that contain (color-superconducting) quark matter can undergo 
transients of this type. This is the main difference of our model 
to the hadronic alternatives which work also for low (or canon¬ 
ical 1.4 Mq) stars. It is tempting to think that the various mod¬ 
els of cooling of CS in Cas A can be distinguished through the 
measurements of the mass or other integral parameters of the 
CS in the future. 

In addition to the Cas A behaviour we also studied the ef¬ 
fect of variations of the mass of the CSs on their cooling be¬ 
haviour. First, we find that larger (close to the maximum) mass 
stars will not show the same transient as the assumed model 
with 1.93 Mq mass; rather they will cool to lower temperatures 
much faster. As our code does not cover the very early evo¬ 
lution of CSs, we cannot resolve the transient that takes place 
very early in the evolution of a highly massive CS. Secondly, 
we find that low mass stars M ~ \.\ Mq are much hotter than 
their heavy counterparts, so that the band covered by the stars 
in the range 1.1 < M/Mq < 2 can account for the measured 
surface temperatures of thermally emitting CS. 

There are a number of issues that need further attention; 
(a) The appearance of strangeness in form of hyperons in the 
hadronic phase and strange quarks in the quark phases can sub¬ 
stantially affect cooling. In particular strange quarks may even 
change the structure of the phase diagram of color supercon¬ 
ducting matter at high density by inducing new phases and new 
channels for neutrino emission see, for example. Ref. ll3l];(b) 
heating in the core and crust of the star can become important at 
later stages of evolution; it does not affect the cooling at the age 
of CS in Cas A (if the CS is not a magnetar); (c) the density de¬ 
pendence of pairing gaps of red-green and blue condensates can 
be a factor. The density of states at the Fermi surface increases 
as 1/3 power of density, whereas the strong coupling constant 
decreases logarithmically. Therefore, one can envision a situ¬ 
ation where the red-green condensate has ^ < 1 only in part 
of the core. This will reduce the effective volume of the quark 
phase contributing to the fast cooling; (d) the superconducting 
quark phases may gradually appear as a result of deceleration 
of the star, i.e., the CS can experience a rotation induced phase 
transition to a color superconducting phase iH . This transition 
will occur in the densest region (center of the CS) in parallel 
with deconfinement of quarks. 
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